Tests of stellar model atmospheres by optical interferometry 

VLTI/VINCI limb-darkening measurements of the M4 giant if/ Phe* 
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Abstract. We present /f-band interferometric measurements of the limb-darkened (LD) intensity profile of the M 4 giant star 
if/ Phoenicis obtained with the Very Large Telescope Interferometer (VLTI) and its commissioning instrument VINCI. High- 
precision squared visibility amplitudes in the second lobe of the visibility function were obtained employing two 8.2 m Unit 
Telescopes (UTs). This succeeded one month after light from UTs was first combined for interferometric fringes. In addition, 
we sampled the visibility function at small spatial frequencies using the 40 cm test siderostats. Our measurement constrains 
the diameter of the star as well as its center-to-limb intensity variation (CLV). We construct a spherical hydrostatic PHOENIX 
model atmosphere based on spectrophotometric data from the literature and confront its CLV prediction with our interferometric 
measurement. We compare as well CLV predictions by plane-parallel hydrostatic PHOENIX, ATLAS 9, and ATLAS 12 models. We 
find that the Rosseland angular diameter as predicted by comparison of the spherical PHOENIX model with spectrophotometry is 
in good agreement with our interferometric diameter measurement. The shape of our measured visibility function in the second 
lobe is consistent with all considered PHOENIX and ATLAS model predictions, and significantly different from uniform disk 
(UD) and fully darkened disk (FDD) models. We derive high-precision fundamental parameters for if/ Phe, namely a Rosseland 
angular diameter of 8.13 ± 0.2 mas, with the Hipparcos parallax corresponding to a Rosseland linear radius R of 86 ± 3 R Q , 
and an effective temperature of 3550 ± 50 K, with R corresponding to a luminosity of logL/L o =3.02 ± 0.06. Together with 
evolutionary models, these values are consistent with a mass of 1.3 ± 0.2 M Q , and a surface gravity of logg = 0.68 ± 0.11. 
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1. Introduction 

Stellar atmosphere models predict the spectrum emerging from 
every point of a stellar disk. However, model atmospheres are 
usually only constrained by comparison to integrated stellar 
spectra. Optical interferometry has proven its capability to go 
beyond this principal test of the predicted flux, and to probe the 
wavelength-dependent center-to-limb intensity variation (CLV) 
across the stellar disk. In addition, the measurement of the stel- 
lar angular diameter together with the bolometric flux is the pri- 
mary measure of the effective temperature, one of the most im- 



to: 



M. Wittkowski, e-mail: 



Send offprint requests 
mwittkow@eso . org 

* Based on public data released from the European Southern 
Observatory VLTI obtained from the ESO/ST-ECF Science Archive 
Facility. The VLTI was operated with the commissioning instrument 
VINCI and the MONA beam combiner. 



portant parameters for modeling stellar atmospheres and stellar 
evolution. 

Further tests of stellar atmosphere models by interferomet- 
ric observations help to improve the reliability of results in all 
areas of astrophysics where such models are used. Cool giants 
and AGB stars are of interest for atmosphere modeling since 
they allow the study of extended stellar atmospheres and the 
stellar mass-loss process. Red giants are also used as probes of 
the chemical enrichment history of nearby galaxies through de- 
tailed abundance measurements of the Calcium infrared triplet 
which rely on model atmospheres. 

However, the required direct measurements of stellar in- 
tensity profiles are among the most challenging programs in 
modern optical interferometry. Since more than one resolu- 
tion element across the stellar disk is needed to determine 
surface structure parameters beyond diameters, the long base- 
lines needed to obtain this resolution also produce very low 
visibility amplitudes corresponding to vanishing fringe con- 



trasts. Consequently, direct interferometric limb-darkening ob- 
servations of stars with compact atmospheres, i.e. visibility 
measurements in the 2nd lobe, have so far been limited to a 
small number of stars (including Hanbury Brown et al. 119741 
Di Benedetto & Foy [19861 Quirrenbach et al. [19961 Burns 
et al. P1W71 Hajian et al. ITWgl Wittkowski et al. EUoTY For 
stars with extended atmospheres, described by for instance 
Gaussian-type or two-component-type CLVs, measurements 
of high to medium spatial frequencies of the 1st lobe of the 
visibility function may lead to CLV constraints as well (see, 
e.g. Haniff et al. 1995 Perrin et al. 1999 1. Furthermore, lu- 
nar occultation measurements may also enable a reconstruc- 
tion of the CLV (e.g. Richichi et al. I1990> . Recent optical 
multi-wavelength measurements of the cool giants y Sge and 
BY Boo (Wittkowski et al. 2001 1 succeeded not only in di- 
rectly detecting the limb-darkening effect, but also in con- 
straining ATLAS 9 (Kurucz 1993 1 model atmosphere parame- 
ters. Aufdenberg & Hauschildt ( 2003 1 showed that these y Sge 
interferometric measurements and ySge spectroscopic mea- 
surements both compare well with predictions by the same 
spherical PHOENIX model atmosphere (Hauschildt et al.[l999 1. 

Here, we present limb-darkening observations of the M4 
giant star t/^Phe (HD 11695, HR555, FK5 67, HIP 8837), ob- 
tained during the commissioning period of the ESO Very Large 
Telescope Interferometer (VLTI) with its commissioning in- 
strument VINCI. VINCI is operated with one near-infrared 
/f-band filter and could not provide measurements at differ- 
ent spectral bands. Spectrally resolved VLTI measurements are 
planned with the upcoming instruments AMBER (Petrov et al. 
2003 1 for near-infrared wavelengths, and MIDI (Leinert et al. 
2003 1 for mid-infrared wavelengths. We construct a spherical 
hydrostatic PHOENIX model atmosphere based on spectropho- 
tometric data from the literature. This is the usual procedure 
to obtain a model atmosphere since additional interferomet- 
ric observations are usually not available. This resulting at- 
mosphere model predicts the LD intensity profile, a predic- 
tion that we confront with our interferometric measurement 
in order to test it. Agreement of the model prediction and our 
measurement increases confidence in atmosphere modeling for 
cool giants. In addition, we compare predictions of a plane- 
parallel PHOENIX model atmosphere as well as of standard 
plane-parallel ATLAS 9 (Kurucz [19931 ) and ATLAS 12 (Kurucz 
1996l l2003> atmosphere models. 

The observations and methods presented here are also a 
precursor for more detailed limb-darkening measurements and 
observations of other stellar surface features. The feasibility 
to derive stellar surface structure parameters beyond diameters 
and limb-darkening using VLTI and its scientific instruments 
has been studied by, e.g., von der Liihe ( 1997 1, Jankov et al. 
( I200T1 . and Wittkowski et al. ( 120021 . 

2. Characteristics of </> Phe 

2.1. if/ Phe - A spectroscopic binary ? 

The Hipparcos Catalogue (Perryman & ESA 1997 1 notes in 
the spectral type column that i/'Phe is a spectroscopic binary 
("M4III SB"). (A Phe is hardly worthy of the spectroscopic 



binary designation. Ihere is no mention of this star in the 
"Catalogue of the orbital elements of spectroscopic binary sys- 
tems" (Batten et al. 119891. Likely sources for the "SB" note 
are Houk d!978i and Wilson (I1953I I. The radial velocity from 
Wilson ( 1953 1 is noted with a (poor) grade "c" and the desig- 
nation "SB". The "SB" appears to be based on a note by Lunt 
dl919b : "either obviously variable in velocity or suspected". 
Since then, a small number of radial velocity measurements 
(Evans et al.l!957l Jones fT972l Crampton & Evans ll973l Jones 
& Fischer 1984 1 have failed to show variations in if/ Phe's radial 
velocity. 

While the variability of if/ Phe's radial velocity has not 
been convincingly established, it is clear that if/ Phe is a small- 
amplitude (AV = +0.1), short-period (^ 30 d), photometric 
variable star (Eggen [T973l ). Our high spatial frequency obser- 
vations (see Sect. 14. II cover two days or ~ 1/15 of the pho- 
tometric period. We are not able to establish the photometric 
phase of the interferometric observations since there was no 
follow-up to Eggen's original photometry. 

2.2. The bolometric flux of ip Phe 

We determine the bolometric flux of i/'Phe to be Fb i = 
(3.2 + 0.3) x 10~ 9 WirT 2 , derived from a spline fit and in- 
tegration of the available spectrophotometry. We have used 
the narrow-band spectrophotometric data (A/l = 2.5 nm, 400 
nm - 700 nm) from Burnashev ( 1985) and the medium- and 
broad-band optical/near-IR photometric data from Johnson & 
Mitchell 119751 WTh and Feast et al. dl990l . The absolute 
calibration of the 13 -color photometry is taken from Johnson 
& Mitchell JT9751 [19771 . T he abs olute calibration of JKLM 
photometry is from Johnson J1965b while the H-band calibra- 
tion is from Bessel & Brett ( 1988). We have assumed a uni- 
form 5% error for all absolute fluxes. The color excess towards 
i/'Phe is very low due its high galactic latitude (b = -67°). The 
COBE dust maps (Schlegel et al. 1998) indicate that the maxi- 
mum E(B - V) value for this line-of-sight is 0.026 and we take 
this as the upper limit on E(B - V). This uncertainty in the in- 
terstellar extinction is included in the derived bolometric flux. 
The V-band variability of if/ Phe, which we have not included, 
may be a source of additional uncertainty in F\, \. 

2.3. A priori estimate of fundamental parameters 

In order to decide on the grid of if/ Phe's mass, radius, and ef- 
fective temperature to explore, we roughly estimate these pa- 
rameters. The diameter of i/'Phe has never been measured be- 
fore by a direct technique. Using the calibration of Dyck et 
al. (1996), based on the spectral type and the -magnitude, 
we derive a Rosseland angular diameter of 8.0 + 0.8 mas. 
Here, the Rosseland radius is defined as the radius at which 
the Rosseland optical depth equals unity, a definition that we 
follow in this article. 

With our value for F\, \, T e ff is constrained to 3500 + 260 
K. With the Hipparcos parallax, the linear Rosseland radius is 
derived to 85 + 10 R Q , and the luminosity to log LfL Q - 2.99 + 
0.23. With these values we can place i/'Phe on the theoretical 



lable 1. Properties or i/'Fhe. Ihe upper table lists the adopted literature values (visual magnitude, near-infrared A-band mag- 
nitude, spectral type, metallicity, parallax, bolometric flux). The lower table compares the Rosseland angular diameter and the 
effective temperature with their 1 cr errors as derived in this work using different methods, as well as derived parameters (effective 
temperature, Rosseland radius, luminosity, mass, surface gravity) with their errors. 



K 



Sp. Type [Z/Z ] 



n 



bol 



i/rPhe 4.3-4.5 (period ^ 30 d) 1 -0.63 ± 0.02 2 M4 III 3 0.25 ±0.1 4 
1 Eggen (T573); 2 Gezari et al. (19991 : 3 Houk (19781 ; 4 Feast et al. (19901 



10.15 ±0.15 mas 5 (3.2 ± 0.3) x 10-" WrrT 2 6 
; 5 Perryman & ESA (19971 : 6 see Sect.l2~2l 
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estimates 
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Fit of 

spherical spherical 

PHOENIX model PHOENIX 
to spectrophotometry 

Sect.fJH Sect.|5l 



Fit of 

plane-parallel plane-parallel plane-parallel 

PHOENIX ATLAS 9 ATLAS 12 

model to our interferometric data 

SectH SectH SectJU 



Final 
values 



Sect. 



Ross (mas) 
T eB (K) 

r cff (K) 

R/Re 
log LI L G 
M/M Q 
logg (cgs) 



8.0 ±0.8" 



3500 ± 260 h 

85 ± 10 c 
2.99 ± 0.23 d 

i 9 +0.8 e 
1 - Z -0.6 

66 +0 ' 33 f 

uuu -0.57 



8.0 ±0.4 
3550 ± 50 



8.13 ±0.2 



8.17 ±0.2 



8.24 ± 0.2 



8.19 ±0.2 



8.13 ± 0.2 s 
3550 ± 50 * 
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: 0.08 d 
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3472 
86 

3.02: 



t 125 " 

: 0.06 d 



66 +01S / 
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1.3 ±0.2'' 
0.68 W 



° Using the calibration by Dyck et al. 1 1996) based on the spectral type and the A*-magnitude; b With #r oss and F\, \ \ c with #r oss and n; d with R 
and r cff ; ' with L, r cff , and the evolutionary tracks by Girardi et al. ( 2000 1, see Fig.Q f with M and R, s from spherical PHOENIX model fit to 

interferometric data, from spherical PHOENIX model fit to spectrophotometry. 



H-R diagram (see Fig^. Together with the evolutionary tracks 
of Girardi et al. (J2000 1, these values for r e fr and L give a mass 
estimate of 1.2 ^ M Q , and a logg estimate of 0.66 t^'sr Feast 
et al. ( fT990b derive a metallicity [Z/Z ] of 0.25 + 0.08, based 
on JHK photometry and atmosphere models by Bessell et al. 
( 1989 1. However, they argue that this offset might be an artifact 
since [Z/Z ] = should be appropriate for local M giants. The 
values given in this section are summarized in Tabled 

3. Calculation of atmosphere models for t/y Phe 

The modeling of atmospheres of cool giant stars, such as 
if/ Phe, is complicated by two effects, the treatment of molec- 
ular opacities (molecules form due to the low temperatures), 
and the effects of their spherical extension (cf. Plez et al. 1992 
Hauschildt et al. 1999 1. The appropriateness of using a plane- 
parallel model is based on the photospheric scale height rela- 
tive to the radius of the star. For the sun, with a photospheric 
thickness of ~ 1000 km, the extension of the photosphere rel- 
ative to the solar radius is ~0.1%. For i/rPhe the predicted 
photospheric thickness is ~ 4,000,000 km and with a radius of 
~ 857?0, the extension is ~6%. Atmosphere models for M gi- 
ants including spherical extension effects have been discussed 
by, e.g., Scholz & Tsuji JT984b ; Scholz ( fl985l ; Bessell et al. 
J1989I199U : Plez etal. fl992l : Hofmann & Scholz (fT998l : and 
Hauschildt et al. (1999 1. Here, we employ hydrostatic spher- 
ically symmetric PHOENIX models including state-of-the-art 
treatment of molecular opacities. These models are based on 
the PHOENIX/NextGen models by Hauschildt et al. dl999> . 
We compare our observational results as well to a hydrostatic 
plane-parallel PHOENIX model. Tabulated model CLVs based 
on hydrostatic plane-parallel ATLAS 9 models are publicly 



available (Kurucz 1993 1 and often used for the interpretation of 
interferometric measurements of different stars including cool 
giants. These models are also used to parametrize LD profiles 
(for instance, Claret l2000ll2003l Davis et al. l2000b . We com- 
pare the predictions by these standard models as well. Finally, 
we use a hydrostatic plane-parallel ATLAS 12 model (Kurucz 
[T996ll2003l which, compared to the ATLAS 9 models, includes 
an improved treatment of the molecular opacities, but no spher- 
ical extension effects. 

3.1. PHOENIX model atmospheres 

We have computed new, fully line-blanketed (7xl0 5 atomic 
lines and 9xl0 7 molecular lines), spherical, hydrostatic atmo- 
sphere models with solar photospheric abundances (Grevesse 
& Noels fl993l with version 13 of the PHOENIX code (for a 
general description see Hauschildt & Baron 1999 1. The three 
most important input parameters for our spherical models are 
the effective temperature, the surface gravity, and the stellar 
mass. Our new models start from the NextGen model grid of 
Hauschildt et al. ( 1999 1 with a stellar mass of 0.5 M . The mi- 
crotubulence for all our new models is 2kms~', as adopted for 
the NextGen grid. The new model structures are iteratively con- 
verged with temperature corrections following energy conser- 
vation for convective and radiative equilibrium. We have com- 
puted a grid of 63 new models at seven temperatures: 3450, 
3475, 3500, 3525, 3550, 3575, 3600 K; three gravities: 0.5, 
0.7, 1.0; and three masses 0.7, 1.0, 1.3 M . In addition we have 
computed 27 models at masses 2.0, 2.5, 3.0 M for temper- 
atures 3500, 3550, 3600 K and gravities 0.5, 0.7. 1.0 in or- 
der to study the effect of higher masses; as well as 7 addi- 
tional models at gravities 0.0, 0.25, 0.5, 0.7, 1.0, 1.25, 1.5 for 



mass 1.0 M and temperature 3550 K in order to study the ef- 
fects of different gravities over a larger range. For each model, 
we tabulate the intensity profile at 64 angles for wavelengths 
from 1.8yum to 2.5 jum in steps of 0.5 nm. The mass and in- 
put surface gravity fix the reference radius (at r(Rosseland)=l) 
at which the effective temperature is defined. Furthermore, the 
mass effectively regulates the relative extension of the atmo- 
sphere for fixed values of the effective temperature and grav- 
ity. In hydrostatic equilibrium, the thickness of the atmosphere 
relative to the overall radius increases as the mass decreases. 
One way to quantify the relative extension or degree of com- 
pactness of a hydrostatic stellar atmosphere is to compare the 
gas-pressure scale height with the stellar radius (see Bessell et 
al. 119911 Baschek et al. I199U . In extended atmospheres the 
pressure scale height is a significant fraction of the stellar ra- 
dius, thus the ratio of these quantities is a useful quantity. For a 
fixed gravity, the stellar mass does not directly affect the pres- 
sure scale height, Hp = where % is the gas constant, T 
is the temperature, // is the mean molecular weight, and g e s 
is the effective gravity which includes radiative and turbulent 
acceleration terms. On the other hand, the stellar mass does di- 
rectly affect the stellar radius, R = sjGM/g, where G is the 
gravitational constant, M is the stellar mass, and g is the grav- 
ity. Hence the ratio Hp/R decreases with increasing mass for 
a fixed value of gravity. Therefore, in the hydrostatic context, 
spherical atmospheres for massive stars are more compact than 
for lower mass stars with the same effective temperature and 
gravity. As a result the mass can affect synthetic limb inten- 
sity profiles and therefore can be an important parameter when 
comparing spherical model atmosphere predictions to interfer- 
ometric data. Angular diameter, temperature, gravity, and mass 
are, of course, dependent parameters and can not be indepen- 
dently explored over large ranges. 

For comparison, we use as well a hydrostatic plane-parallel 
PHOENIX model with an effective temperature of 3550 K and 
logg = 0.7. The intensity profile is tabulated at 64 angles for 
wavelengths from 1.8yum to 2.5 //m in steps of 0.5 nm, as for 
the spherical models. 

3.2. Comparison of PHOENIX models with 
spectrophotometry 

The effective temperature and angular diameter of i/fPhe can 
be constrained by comparing synthetic spectral energy dis- 
tributions (SEDs) from the spherical PHOENIX models with 
the observed spectrophotometry. To do so, synthetic photom- 
etry is performed using the 13-color filter sensitivity functions 
from Johnson & Mitchell ( 1975 ) compiled by Moro & Munari 
( 2000 1 and the JHKLM filter sensitivity functions from Bessell 
& Brett (1988 1. The synthetic spectrophotometry values are 
scaled to the observed values by the factor 

^ ''Ross lu 

where 0r oss is the wavelength-independent Rosseland angular 
diameter, A^/Ay is the extinction curve relative to the V-band, 
E(B - V) is the color excess, and we assume Ry - 3.1 = 
Ay I E(B-V) using the extinction curves of Cardelli et al. ( 1989 1 
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Fig. 1. The predicted regions for i/rPhe in the Hertzsprung- 
Russell diagram from the model atmosphere fitting to the ob- 
served spectral energy distribution. The method provides val- 
ues for T e ff and 0r oss , and together with the parallax, a luminos- 
ity. The error ellipses show 1- and 2-sigma zones. Also shown 
are the evolutionary tracks of Girardi et al. (2000 1 for masses 
M > 0.6 M Q . 



A comparison of the best fit synthetic SED to the 13-color 
photometry, which covers the widest wavelength range con- 
temporaneously, with all the spectrophotometric data is shown 
in Fig. 13 The shape of the near-infrared continuum and the 
strong TiO bands in the optical provide the principal constraints 
on T e ff, while the fluxes in absolute units constrain #r oss . The 
best fit parameters, based on least-squares fits to the photom- 
etry and application of the F-test (Aufdenberg et al. 2002) are 
T eS = 3550+ 50 K and 6 Ross = 8.0+ 0.4 mas. Together with the 
parallax (Tabled, the ranges in these values translate into the 
1- and 2-sigma uncertainty regions on the theoretical HR dia- 
gram in Fig. [2 The uncertainty in the luminosity is somewhat 
of an overestimate, since the best fit values for T e ff and 6Ir oss are 
not independent. 

While T e g and 6Ir oss are well constrained by the SED fit 
alone, the surface gravity and mass are not. However, relying 
on the evolutionary tracks from Girardi ( 2000 1, the 2-sigma re- 
gion on the HR diagram indicates that the mass is 1.2 ± 0.4 
M . A surface gravity then follows from this mass, #r oss , and 
the parallax: log g — 0.66 ^ 23- Since these constraints are much 
tighter than can be derived from comparisons to spectropho- 
tometry alone, the mass and surface gravity parameters for the 
synthetic spectrum shown in Fig.[2]were chosen to be consis- 
tent with these values. Table[Ogives a summary of the derived 
values. 

3.3. ATLAS model atmospheres 

In addition to the PHOENIX models, we employ intensity pro- 
files predicted by standard plane-parallel hydrostatic ATLAS 9 
model atmospheres from the Kurucz CD-ROMS (Kurucz 
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Fig. 2. (Top) Spectral energy distribution of iJ/Phe (error bars) compared with synthetic photometry ('x's) for the narrow and 
broad-band photometry derived from a high-resolution model spectrum (spherical PHOENIX model) shown in gray. (Bottom) The 
percentage difference between the observed and synthetic photometry in each wavelength bin. The model parameters are derived 
by a least-squares fit to only the 13 -color photometry (Johnson & Mitchell ll975> . while this Fig. shows the model comparision 
to all available data, as described in the text. 



I1993> . The Kurucz CDROMs tabulate monochromatic LD pro- 
files for 17 angles in 1221 frequency intervals ranging from 
9.09 nm to 160 000 nm. In the range of the near-infrared K- 
band filter (1 800-2500 nm), the frequencies are sampled in 
steps of 10 nm. This data is available for different chemi- 
cal abundances, microturbulent velocities, effective tempera- 
tures and surface gravities. To be consistent, we have chosen 
the parameters to be closest to those determined for our fa- 
vorite PHOENIX model, i.e. the model that best fits the spec- 
trophotometry (Sect. 13 .21 . These parameter values are T e s 
3500 and 3750 K, logg 0.5 and 1.0, solar chemical abundance, 
and a standard microturbulent velocity of 2kmsec~' (file ku- 
rucz. cfa.harvard.edu/grids/gridP00/ip00k2.pckl9). 

Finally, we employ a plane-parallel hydrostatic ATLAS 12 
model which includes an improved treatment of the opacities 
(Kurucz 1996, 2003 1. We have chosen again the parameters of 
our favorite spherical PHOENIX model for reasons of consis- 
tency, i.e. T e ff =3550 K, logg=0.7, and solar chemical abun- 



dance. The ATLAS 12 intensity profile is tabulated at 16 angles 
for wavelengths from 1.8 fim to 2.5 fim in steps of 0.5 nm. 

3.4. Effects of model geometry on the interpretation of 
fitted angular diameters 

We use spherical PHOENIX models as well as plane-parallel 
PHOENIX and ATLAS models, as described above, for the fits 
to our interferometric data. The fit result for both geometries 
is the LD diameter 6*ld at which the intensity reaches zero 
(see Sect. 14.51 below). Figure [3] illustrates CLV predictions 
by spherical and plane-parallel PHOENIX models (7^=3550 K, 
logg=0.7, mass of the spherical model 1.3M ; CLVs are cal- 
culated for the broad VINCI K passband). The monochromatic 
contributions over the /T-band to this passband-averaged spher- 
ical CLV are shown below in Fig. [5] (middle). It illustrates 
that our broad-band measurement is dominated by continuum 
/T-band photons, and that it has only little line contamination 
(relative to a strong molecular band). The escaping continuum 
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Fig. 3. Comparison of CLVs predicted by a spherical (solid 
lines) and a plane-parallel (dashed lines) PHOENIX model, both 
with the same 0% intensity diameter #ld- Both models have 
parameters T e ff =3550 K, logg=0.7. The spherical model has a 
mass of 1.3 M Q . Both CLVs are computed for the broad-band 
VINCI sensitivity function (A'-band). The two CLVs differ by 
a tail-like extension for the spherical model which is missing 
for the plane-parallel model. We describe the diameters of both 
models by the Rosseland diameter, which is close to the point 
where the CLV drops steeply and most continuum photons es- 
cape. The f?R OSS values for these two model geometries differ by 
a factor of 0.94, while the 9 LD values correspond. 



photons originate from a rather compact zone at /j ~ 0.3. The 
spherical CLV shows an inflection point and a tail-like exten- 
sion which is caused by the optically-thin limb of the spherical 
model. The plane-parallel model is semi-infinite for all angles 
6 and has a singularity at 6 — 90° or = 0m. This causes the 
sharp edge of the plane-parallel CLV and its lack of any tail- 
like extension. As a result, despite the correspondence of the 
#ld values, the points where the CLVs drop steeply, which is 
also the point where most continuum photons escape, do not 
correspond for these two model geometries. In other words, dif- 
ferent 6*ld fit values are expected for plane-parallel and spheri- 
cal model geometries using the same physical data. We use the 
Rosseland diameter c9r oss , which is associated with a Rosseland 
optical depth of unity, to describe the stellar diameters. The 
Rosseland diameter is close to the compact zone where most 
continuum photons escape and where the CLV drops steeply, 
and is hence comparable for both model geometries. 
For the plane-parallel models, we assume 0m = #Ross, since the 
plane-parallel CLVs drop steeply to zero directly at the stellar 
limb. 

For the spherical models, 0m is associated with the outermost 
shell of the model (at Rq) and the wavelength-independent con- 
version factor from 0m to #r oss is the ratio Cr oss /ld = R(tr ss = 
1)/Rq. Here, the outermost radius Rq of the model is defined 
by the standard outer boundary conditions for our model at- 
mospheres, which are a continuum optical depth of le-6 at 
1.2yum and an outer gas pressure of le-4 dynes/cm 2 . The in- 
ner optical depth boundary is 100 at 1.2 /mi. The radial grid 
and outermost radius differ from model to model, resulting 
in a model-dependent conversion factor Cr oss /ld (as used in 
Table|U Sect.|51below). 



4. vlii/vinui measurements 

4.1. Observations 

The ESO VLTI, located on Cerro Paranal in northern Chile, is 
currently in the phase of commissioning. For a recent general 
overview see Glindemann et al. ( 2003 ) and references therein. 
First fringes using 40 cm test siderostats and the commission- 
ing instrument VINCI (Kervella et al. 2003a I were achieved in 
March 2001, followed by first fringes using the 8.2 m diame- 
ter Unit Telescopes (UTs) on October 30, 2001. First fringes 
with the first scientific instrument, the mid-infrared instrument 
MIDI (Leinert et al. 2003 1, were obtained in December 2002. 

Our observations of i^Phe's squared visibility amplitudes 
beyond the first minimum employing the UTs succeeded one 
month after the achievement of first fringes with the UTs, and 
were first mentioned in ESO press release 23/01 (Nov. 5, 2001). 
These visibility values of t// Phe were obtained during the nights 
starting on Nov. 1 & 2, 2001 using the UT 1 - UT 3 102 m base- 
line of the VLTI, and the VINCI instrument, as part of the VLTI 
commissioning. The UTs were not yet equipped with adaptive 
optics nor a tip-tilt corrector. High-contrast fringes at short spa- 
tial frequencies were obtained before and after the UT measure- 
ments using the 40 cm test siderostats. These data, as all scien- 
tifically interesting VLTI/VINCI data, have been made publicly 
available through the ESO archive. Table |2] details our obser- 
vations and Table [3] shows the adopted properties of the used 
calibration stars. 

The use of the large apertures of the 8 m UTs together 
with the use of single-mode optical fibers, a technique that is 
known to lead to very accurate visibility measurements (see, 
e.g. Coude du Foresto 1998), has enabled precise measure- 
ments of low visibility amplitudes. Bootstrapping techniques, 
such as those used by Hajian et al. {1998 ) and Wittkowski et 
al. (|2001 ) to detect and track fringes with weak contrasts, were 
not necessary. 

Configurations using different aperture sizes for different 
baseline lengths, as used here, might be an interesting per- 
spective for the completed VLTI as well. The completed VLTI 
will routinely allow the measurement of the limb-darkening 
effect for many stars, and open the possibility to obtain stel- 
lar surface structure parameters beyond limb-darkening (see 
e.g., Wittkowski et al. 2002 1. This will require visibility data 
uniformly distributed over the wv-plane including both, high- 
contrast fringes in the first lobe as well as very low-contrast 
fringes in the second and following lobes of the visibility func- 
tion. Then, the lowest fringe contrasts can be measured with 
the 8 m UTs in order to obtain a sufficient signal-to-noise ra- 
tio (S/N) for fringe detection and tracking and/or an enhanced 
S/N for these important low visibility points while higher fringe 
contrasts on the same source can be obtained with the 1.8 m 
ATs. 

Our data were recorded through the VINCI A"-band filter. 
The scan length of the optical path difference was 280 fan, the 
fringe frequency, i.e. the time needed to scan one interferomet- 
ric fringe, was 296 Hz for the siderostat data and 695 Hz for the 
UT data. Our data were taken as series of 100 or 500 interfero- 
gram scans. 



lable 2. Record or our t/rFhe observations. Listed are the date 
of observation, the used stations, the baseline length, the num- 
ber of series of interferograms, and the used calibration stars. 



Date 




Stations 


S(m) 


# 


Calibration stars 


2001- 


10-12 


E0/G0 


16 


4 


ySCet, eLep 


2001- 


10-16 


E0/G0 


16 


3 


ySCet, r/Cet 


2001- 


10-17 


E0/G0 


16 


3 


jSCet, 7/Cet 


2001- 


10-18 


E0/G0 


16 


3 


/JCet, ^Cet 


2001- 


11-01 


UT 1/UT 3 


102 


13 


^•Phe, y 2 Vol 


2001- 


11-02 


UT 1/UT 3 


102 


5 


39Eri, y 2 Vol 


2001- 


11-04 


E0/G0 


16 


7 


/JCet, eLep 


2001- 


11-05 


E0/G0 


16 


3 


/JCet, eLep 



Table 3. Properties of the used calibration stars. Listed are the 
spectral type from Perryman & ESA ( 1997 1, the /f-band magni- 
tude from Gezari ( 1999 1, the effective temperature T e ff derived 
from the spectral type, the effective wavelength Aq derived from 
r e ff and the VINCI sensitivity curve, the adopted /f-band UD 
diameter and its error eQ. The sources for the adopted diameter 
values and their errors are (1) Cohen et al. (1999 1 and/or (2) 
Borde et al. J20021. 



Obj. 


SpT 


K 


Teff 

K 


/im 


mas 


e6 
mas 


Ref. 


/3Cet 


K0III 


-0.2 


4800 


2.179 


5.18 


0.06 


1,2 


eLep 


K4III 


-0.2 


4075 


2.181 


5.91 


0.06 


1.2 


T^Cet 


K2III 


0.9 


4600 


2.180 


3.35 


0.04 


1.2 


^Phe 


K5III 


1.5 


4000 


2.181 


2.69 


0.03 


1.2 


y 2 Vol 


K0III 


1.5" 


4800 


2.179 


2.44 


0.06 


1 


39Eri 


K3III 


2.3 


4200 


2.181 


1.81 


0.02 


1 



from V and the V - K color for a K0 III star. 



4.2. Data reduction 

We have used a modified version (Kervella et al. 2003b I of the 
standard VINCI data reduction pipeline (Ballester et al. l2002> . 
whose general principle is based on the original algorithm of 
the FLUOR instrument (Coude du Foresto et al. 119971 Coude 
du Foresto et al. 1998 1. This data reduction software includes 
the following procedures. 

For each scan, the two interferograms produced from the 
two interferometric outputs of the VINCI beam combiner are 
calibrated for photometric intensity fluctuations using the two 
photometric outputs. After this first photometric calibration, 
the two calibrated interferograms are subtracted to remove 
residual photometric fluctuations. Since the two fringe pat- 
terns are in perfect phase opposition, this subtraction removes 
a large part of the correlated fluctuations while enhancing the 
interferometric fringes. This is particularly important for the 
small tfr Phe visibility amplitudes in the second lobe of the vis- 
ibility function obtained with the UT1-UT3 baseline. Instead 
of the classical Fourier analysis, the time-frequency analysis 
(Segransan et al. 1999 1 based on the continuous wavelet trans- 
form (Farge 1992 1 is used to derive the wavelet power spectral 
density (wavelet PSD) of each scan. The differential piston cor- 
rupts the amplitude and the shape of the fringe peak in the PSD. 
A selection of single scans is used to remove interferograms 
that are affected by strong differential piston and, if existent, 
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Fig. 4. Average wavelets PSD of one series (Date 2001-11-02, 
Univ. Time 05:23:23) of 414 accepted interferograms obtained 
on i^Phe (UT1-UT3 baseline). The thin curve is the simple av- 
erage of all PSD before removal of the background. The back- 
ground is estimated by a fit of a linear function to low and 
high frequency PSD values outside the fringe signal (ranges 
[HOO.^OOOlcirT 1 and [SOOO.^SSJcirr 1 ) and compensated. 
The squared coherence factor is derived from the integration 
of the background calibrated PSD (thick line) in the range 
[2000. . 8000] cm -1 . Before averaging, each fringe peak is re- 
centered to the theoretical frequency of the fringes in order to 
reduce the energy spreading due to the differential piston ef- 
fect before the power peak integration. This is the reason why 
the shape of the non-re-centered (thin line) and re-centered and 
bias corrected (thick line) PSDs are not exactly the same. 

that show no fringes. The selection criteria include checks of 
the peak width in the time domain (+50% around theoretical 
value is accepted), as well as the peak position (+30% around 
theoretical value is accepted) and peak width in the frequency 
domain (±40% around theoretical value is accepted). In addi- 
tion, fringe scans with low photometric signal (S/N of less than 
5), with a large jump of the optical path difference before the 
scan (larger than 20yum), and with fringes at the edge of the 
scan are rejected. The average wavelet PSD of the selected in- 
terferograms is computed. 

The use of squared quantities, such as the PSD, requires at- 
tention to a possible additive background of the PSD caused by 
residual photon and detection noise, despite the first compen- 
sations mentioned above. The relative contribution of remain- 
ing biases would be highest for small visibility values, such as 
those studied here (cf. Wittkowski et al. 2001 1. In the case of 
the VINCI data reduction, the residual noise background is de- 
termined from a least squares fit to low and high frequency PSD 
values outside the fringe signal, and compensated. Figure |4] 
shows an example of the average wavelet PSD derived from 
one series of 414 accepted interferograms obtained on i^Phe 
using UT1 and UT3. Shown are the raw PSD before compensa- 
tion of the residual PSD background, the derived model of the 
PSD background, and the compensated PSD. Despite the very 



strong photometric Fluctuations that are observed in the multi- 
speckle regime, and despite the very low fringe signal (for this 
particular PSD, the squared visibility amplitude of if/ Phe was 
only \V\ 2 = 0.013), the background calibrated wavelet PSD is 
free of any photometric contamination. This is very important 
to reliably obtain our small if/ Phe squared visibility values be- 
yond the first minimum. 

The coherence factor values are finally derived by integrat- 
ing the compensated average wavelet PSD of the selected inter- 
ferograms. This integration includes fringe power at all wave- 
lengths over the broad VINCI sensitivity band. Single scans 
with coherence factors that differ by more than 3 <x from the 
median of the sample are rejected. 

4.3. Calibration 

Table|2]lists the calibration stars that have been used during the 
different nights of our (/'Phe observations; Table[3]details their 
characteristics including the adopted diameters and their errors. 

Owing to the low limiting magnitude of the siderostats 
(K CO n - 1 at the time of these observations), the calibra- 
tion stars for the 16 m siderostat baseline have diameters, be- 
tween 3 mas and 6 mas, which are not small compared to the 
expected if/ Phe diameter of ^8 mas. Thus, the calibration re- 
quires special attention to the adopted diameter values. Cohen 
et al. ( 1999 1 have used a spectro-photometric calibration to de- 
rive high-precision Rosseland diameters and their errors for a 
list of 422 sources. They found that their estimates compare 
well with several other diameter estimates. This was confirmed 
by Borde et al. (2002i, who in addition reduced this list to 
374 stars carefully selected to be used as calibration stars for 
long baseline stellar interferometry. We rely on these diame- 
ter estimates and their errors as given in the lists by Cohen 
et al. (119991 and Borde et al. (20021. The main calibration 
star for the siderostat data, BCet, used during all siderostat 
nights, was in August 2002 calibrated against the Cohen stars 
20Cet (M0 III, 6^ D = 3.41 + 0.037mas) and tCet (K1III, 
6^ D = 3.27 + 0.036 mas) using the 64m G1-E0 baseline. The 
measured B Cet diameter was found to be well consistent with 
the Cohen value and its error given in Table An additional 
calibration uncertainty arises from the fact that the 16 m base- 
line measurements, corresponding to squared visibility ampli- 
tudes in the range ~ 0.8-0.95, only marginally resolve our ob- 
ject i/fPhe. To take these calibration uncertainties for the 16 m 
baseline into account, we weight the measurements obtained 
with this baselines by a factor of two lower than the 102 m 
baseline measurements. A VLTI program in order to derive 
calibration star diameters in a self-consistent way from VLTI 
observations is in progress (Percheron et al. 2003 1. 

For the i/fPhe data sets as well as all calibration star data 
sets, the coherence factors c 2 are computed as described in 
Sect. 14.21 The transfer function f 2 is computed for all calibra- 
tion star data sets as t 2 = c 2 /|V| 2 , where \V\ 2 is the adopted 
squared visibility amplitude for the respective time and base- 
line using the adopted diameter value and effective wavelength 
from Table [3] The transfer function value for the time of each 
of the if/ Phe observations is then obtained as a weighted av- 



erage of the calibration stars t values during the night with 
the weight being a product of (1) a Gaussian function of the 
time difference from the respective if/ Phe measurement with a 
width <j of 3 hours for the siderostat data and 2 hours for the UT 
data, and (2) the error of the single transfer function measure- 
ment consisting of the statistical c 2 error and the error resulting 
from the adopted diameter uncertainty. The final iff Phe squared 
visibility values are then obtained by division of the if/ Phe co- 
herence factors by the transfer function values obtained for the 
respective time. The final errors of the if/ Phe squared visibil- 
ity amplitudes are computed from the error of the if/ Phe co- 
herence factor (the scatter of the single scan's coherence fac- 
tors) and the uncertainty of the transfer function values, which 
include the errors of the calibration stars' coherence factors, 
the adopted diameter errors, and the variation of the computed 
transfer function over the night. 

4.4. Results 

Table |4] shows the resulting squared visibility amplitudes \V\ 2 
for each of our data sets. The given errors of \V\ 2 include the 
statistical error obtained from the scatter of the single scan 
values, the uncertainty of the calibration star's diameter (see 
Table [3}, and the variation of the computed transfer function 
over the night (see Sect. 14. 3> . A plot (Figure [6} of our mea- 
sured | V| 2 values together with best fitting PHOENIX and ATLAS 
models is discussed in Sect.|5](below). 

4.5. Computation of broad-band model visibility values 

The use of the VINCI broad-band K filter influences the ob- 
tained squared visibility amplitudes in two ways (cf. Tango 
& Davis 2002 1. Firstly, the stellar intensity profile is not 
observed at a monochromatic wavelength A but integrated 
over the broad-band sensitivity curve. The strength of the 
limb-darkening effect for tf/Phe, is expected to vary over the 
VINCI sensitivity curve. Figure [5] shows the expected varia- 
tion of i^Phe's flux, intensity profile, and strength of the limb- 
darkening (correction factors from LD to UD diameters) over 
the VINCI passband. The VINCI passband-averaged intensity 
profile and strength of the limb-darkening are shown as well. 
The plots show that our A'-band measurements are dominated 
by continuum photons and that there is only little line contami- 
nation (relative to a strong molecular band). Secondly, squared 
visibility amplitudes are averaged over different spatial fre- 
quencies (B/[(Aq - A/l)..(/l() + A/})]) for any point of obser- 
vation with fixed baseline B and a broad-band filter (central 
wavelength Aq, filter width A A) . This latter effect has been dis- 
cussed by Kervella et al. ( 2003b I for the VINCI instrument and 
been referred to as "bandwidth smearing". 

The calculation of the synthetic broad-band squared visibil- 
ity amplitudes has to follow the way the measured squared vis- 
ibility amplitudes are derived from the raw fringe data. For the 
VINCI data reduction (see Sect. l4~2t . first the PSD, a squared 
quantity corresponding to the squared visibility amplitudes, is 
computed from the interferograms, and then the coherence fac- 
tor is derived from the integral of the PSD, which includes 
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Fig. 5. (Top) Model ^Phe K-band spectrum F A (gray), VINCI 
sensitivity function S ,\ (gray; right scale), and effective spec- 
trum F^Sa (black). The dashed vertical line denotes the ef- 
fective wavelength A® = 2.183yum. (Middle) Monochromatic 
model tfr Phe intensity profiles across the /T-band with a resolu- 
tion of 50 nm (gray) and their S ,\ weighted average (calculated 
with 5 nm resolution, black). With VLTI and AMBER, it will 
in principle be possible to directly measure these different in- 
tensity profiles with a resolution of up to ~ 0.2 nm. (Bottom) 
Monochromatic correction factors from LD to UD diameters 
with a resolution of 5 nm (gray) and 50 nm (black) between 
1.8 /im and 2.5 /urn. The horizontal line denotes the weighted 
average of the correction factors (value 0.917). The spherical 
PHOENIX model as used in Fig.|5]was employed for these plots. 
The spectrum shows CO lines at the red end, which have a sig- 
nificant impact on the tfrPhe intensity profiles at these wave- 
lengths. They become optically thin near the limb and hence, 
the limb is less dark, the radius is larger in the bands, and the 
LD-to-UD correction factor is closer to unity. At most wave- 
lengths the limb edge is at fi ~ 0.3, but at a minority of wave- 
lengths the limb is brighter, almost uniform at one wavelength. 
These profiles come from wavelengths in the CO band. These 
CO lines, however, do not have a strong influence on the VINCI 
passband-averaged model intensity profile because of their nar- 
rowness, their lower intensity, and the shape of S ,\. The LD- 
to-UD correction factors show also the general trend that the 
(continuum) intensity profiles become closer to a uniform disk 
with increasing wavelength. 



lable 4. Resulting squared visibility amplitudes \ V\ tor each 
series of interferograms. Given are the date, the Universal 
Time, the spatial frequency B/Aq (1/") for i/^Phe's effective 
wavelength of /io=2.183/im, the squared visibility amplitudes, 
their error cr(|y| 2 ), and the number (#) of accepted scans. 



Date 


Univ. Time 


B/A 


W\ 2 


(T 


# 


2001- 




an 








10-13 


03:40:05 


35.29 


8.471e-01 


4.902e-02 


462 


10-13 


04:24:29 


35.55 


8.580e-01 


5.038e-02 


439 


10-13 


04:59:41 


35.39 


8.084e-01 


4.938e-02 


460 


10-13 


05:41:03 


34.74 


8.306e-01 


5.054e-02 


452 


10-17 


01:56:26 


33.82 


8.805e-01 


2.306e-02 


479 


10-17 


04:48:29 


35.34 


8.721e-01 


2.650e-02 


450 


10-17 


06:13:42 


33.29 


8.739e-01 


2.742e-02 


483 


10-18 


01:54:08 


33.86 


8.666e-01 


2.082e-02 


463 


10-18 


03:12:03 


35.20 


8.409e-01 


2.270e-02 


395 


10-18 


04:35:33 


35.43 


8.626e-01 


2.442e-02 


433 


10-19 


03:02:01 


35.12 


8.694e-01 


1.754e-02 


444 


10-19 


04:09:31 


35.54 


8.668e-01 


1.712e-02 


481 


10-19 


05:25:09 


34.55 


8.625e-01 


1.739e-02 


464 


11-02 


05:17:37 


201.09 


1.196e-02 


7.545e-04 


50 


11-02 


05:20:05 


200.53 


1.318e-02 


7.164e-04 


398 


11-02 


05:23:23 


199.76 


1.317e-02 


7.220e-04 


414 


11-02 


05:25:21 


199.30 


1.325e-02 


7.388e-04 


62 


11-02 


05:27:49 


198.70 


1.318e-02 


7.309e-04 


413 


11-02 


05:31:00 


197.93 


1.317e-02 


7.371e-04 


411 


11-02 


05:34:13 


197.13 


1.317e-02 


7.425e-04 


433 


11-02 


06:20:40 


183.95 


9.144e-03 


7.091e-04 


66 


11-02 


06:32:42 


180.01 


8.605e-03 


7.322e-04 


77 


11-02 


06:35:13 


179.16 


8.391e-03 


7.009e-04 


308 


11-02 


06:38:38 


178.00 


8.000e-03 


6.891e-04 


329 


11-02 


06:42:02 


176.82 


7.847e-03 


6.887e-04 


286 


11-03 


02:50:05 


221.37 


1.266e-02 


1.385e-03 


81 


11-03 


02:52:31 


221.20 


1.197e-02 


1.360e-03 


387 


11-03 


02:55:47 


220.96 


1.180e-02 


1.369e-03 


354 


11-03 


02:58:52 


220.72 


1.221e-02 


1.410e-03 


386 


11-03 


03:02:10 


220.46 


1.223e-02 


1.430e-03 


426 


11-05 


04:56:32 


33.37 


8.622e-01 


2.822e-02 


348 


11-05 


05:06:35 


32.98 


8.767e-01 


2.458e-02 


492 


11-05 


05:11:46 


32.77 


8.804e-01 


2.466e-02 


492 


11-05 


06:13:25 


29.69 


8.832e-01 


2.218e-02 


195 


11-05 


06:25:05 


29.00 


9.048e-01 


2.592e-02 


413 


11-05 


06:35:29 


28.36 


9.159e-01 


2.572e-02 


453 


11-05 


07:53:25 


23.19 


9.643e-01 


2.640e-02 


427 


11-06 


01:41:53 


35.00 


8.640e-01 


2.294e-02 


454 


11-06 


01:46:19 


35.06 


8.598e-01 


2.312e-02 


484 


11-06 


01:50:41 


35.12 


8.612e-01 


2.336e-02 


479 



fringe power from all wavelengths over the sensitivity band. 
Hence, for VINCI, the model squared visibility amplitudes at 
different wavelengths and spatial frequencies have to be inte- 
grated (see Kervella et al. 2003b), and not the complex visibil- 
ities (visibility amplitudes with sign) as implied by the general 
formalism in Tango & Davis ( 2002 1. This means that the broad- 
band squared visibility amplitude values never reach zero in our 
case. 

We derive non-normalized monochromatic LD model visi- 
bility values Vld(T) on a grid of wavelengths (range 1 .8. .2. 5 fim 
in steps of 0.5 nm for our PHOENIX models) by numerical eval- 
uation of the Hankel transform of the intensity profiles tab- 



ulated tor each or these (1400) monochromatic wavelengths A 
(following Davis et al. 2000 and Tango & Davis 2002 see also 
the use in Wittkowski et al. 2001 1 by 



V LD (A) 



Jo 



S a % Jo [n 6ld (B/A) ( 1 - fi 2 ) 1 / 2 ] p d» 



(1) 



Here, I't is used as a function of //, n = cos 6 (fi e [0,1], yu=0 
corresponds to the star's radius at which the intensity reaches 
zero, and fi-l to the star's center); S a is the VINCI instrument's 
sensitivity curve including the transmissions of the atmosphere, 
the optical fibers, the VINCI /f-band filter, and the detector 
quantum efficiency; B is the sky-projected baseline length. For 
reasons of numerical accuracy, it had been necessary to inter- 
polate the tabulated intensity profiles onto a regularly sampled 
grid of // values (we have chosen a cubic-spline interpolation 
onto 100 n values), rather than interpolating the more complex 
argument of the integral in Eq.^ The integral is then computed 
using a 5 -point Newton-Cotes algorithm. We have checked the 
accuracy of the algorithm by comparing results for tabulated 
UD and FDD intensity profiles with their available analytic ex- 
pressions for V. In order to calculate a spatial frequency B/Aq, 
as used in Fig. [6] (below), an effective wavelength Aq is com- 
puted as 
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with Fa being the stellar flux at wavelength A: 
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(3) 



With our favorite spherical PHOENIX model for i^Phe as de- 
rived in Sect. 13.21 Fig.|2l we obtain Ao = 2.183 //m. This value 
characterizes the wavelength of our measurement, but is not 
used for the calculation of the synthetic visibilities since they 
are calculated for the full wavelength range in Eq.^ The func- 
tions S a and Fa, as well as the effective wavelength are shown 
in Fig. 13 (top). As for the processing of the instrument data, 
each monochromatic visibility is not normalized separately. 
The VINCI passband-averaged squared visibility amplitudes 
are finally computed as 
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(4) 



including the proper normalization by the total detected flux. 

We directly use the model atmosphere's tabulated inten- 
sity profiles (the stellar radiation field) in Eq. ^ N° approx- 
imation of the intensity profile by any limb-darkening law is 
used. This ensures two important points: Parametrizations of 
intensity profiles are usually very good approximations for the 
visibility curve before the first minimum, but may lead to devi- 
ations from the original visibility profile beyond the first mini- 
mum, which may be significant with the precision of our mea- 
surement. Furthermore, the direct use of the tabulated model 
stellar radiation field on a fine grid of wavelengths enables us 
to compute the model squared visibility values for exactly our 



instrument s passband, while published limb-darkening coeffi- 
cients are calculated for specific filter curves which do not well 
resemble the VINCI sensitivity curve including e.g. the trans- 
mission of the optical fibers and the detector efficiency. 

As discussed in Sect. 13.41 the LD diameter used in Eq. ^ 
depends on the detailed model structure and is therefore not a 
well suited quantity for reference. We transform #ld into the 
Rosseland angular diameter # Ross with the factor Cr oss /ld de- 
fined in Sect.l3~4l 



5. Comparison of our VLTI/VINCI data to our model 
predictions 

We use the different model atmospheres as described in Sect. [5] 
and, with the formalism described above in Sect. 14.51 calculate 
synthetic squared visibility amplitudes for the spatial frequen- 
cies of our data. For each considered model, we find by a non- 
linear least squares fit the 6*ld value for which the synthetic 
and measured squared visibility amplitudes have the lowest^ 2 
value. The LD angular diameter 6*ld is treated as the only free 
parameter. The best fitting 6>ld value is then transformed into 
the Rosseland angular diameter #r oss using the factors Cr oss /ld 
defined in Sect. EH 

Table [5] lists the fit results for the different spherical 
PHOENIX, and plane-parallel PHOENIX, ATLAS 12, and ATLAS 9 
models. The first three columns specify the model parame- 
ters r e ff , log g, and M, followed by the model-specific factor 
Cross/ld- The next two columns give the best fitting #ld values 
together with the reduced^ 2 values. Finally, the results for 0r oss 
are given in the last column. The formal errors for 6*ld are found 
to be about uniform for all considered models. The standard 
formal error cr(# LD ) for #ld, derived as the ASld corresponding 
to Ax 2 = 1, is -0.05 mas. Deriving the error by means of the 
F-test, which is a more reliable estimate than the^ 2 -test (e.g., 
Bevington & Robinson 1992), we obtain cr(# LD ) ~ 0.12 mas. 
Here, we assume a number of degrees of freedom of 39, i.e. 
the number of visibility points minus one. The real situation is 
more complicated because (1) the errors include statistical er- 
rors as well as systematic errors (which can not easily be sepa- 
rated), (2) some squared visibility value measurements are not 
independent but linked via the same calibration measurements, 
(3) the visibility function is sampled at two different spatial fre- 
quency regions using different telescopes, and (4) obtained xl 
values are larger than unity. We estimate the final error cr(# LD ), 
including calibration uncertainties, to be ~0.2 mas by com- 
paring results obtained by different calibrations and different 
sub-samples of our data. Table[5]shows results for our favorite 
model atmosphere from Sect. 13.21 i.e. the spherical PHOENIX 
model with r eff =3550K, logg=0.7, and M=1.3M . In addi- 
tion, we consider spherical PHOENIX models with varied pa- 
rameters that are still consistent with the estimates in Sect. 13. 21 
Furthermore, we compare plane-parallel PHOENIX, ATLAS 12, 
and ATLAS 9 models with, for reasons of consistency, parame- 
ters closest to our favorite model. 

Figure [6] plots our measured squared visibility values 
together with the model prediction by our favorite spheri- 
cal PHOENIX model with best fitting 6*ld value. Shown are 
also the predictions by corresponding plane-parallel PHOENIX, 
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Fig. 6. Our measured squared visibility amplitudes of t^Phe ('x' symbols with error bars) together with the (solid black line) 
spherical PHOENIX model prediction with model parameters T e ff, log g, and mass as derived from spectrophotometry and model 
evolutionary tracks (Sect. EOl . and best fitting 6>ld value. Shown are also the (dashed-dotted line) plane-parallel PHOENIX model, 
(dotted line) plane-parallel ATLAS 12 model, (dashed line) plane-parallel ATLAS 9 model, all with corresponding model parame- 
ters and best fitting 6*ld- As a reference for the strength of the limb-darkening, the gray lines denote corresponding UD (upper 
line) and FDD (lower line) model visibility functions. The left panel shows the full range of the visibility function while the right 
panel is an enlargement of the low squared visibility amplitudes in the second lobe. All considered PHOENIX and ATLAS model 
predictions result in a very similar shape of the visibility function in the 2nd lobe. Our measurements are significantly different 
from uniform disk and fully-darkened disk models, and consistent with all considered PHOENIX and ATLAS models. 



Table 5. Results for #ld obtained by fits to our interferometric 
data. Listed are the model input parameters T e ff , log g, and mass 
M, the model-specific correction factors Cr oss /ld, the fit results 
for 6>ld, their corresponding reduced xl values, and finally the 
Rosseland angular diameter 0r oss - The error ct^ld) is estimated 
to be uniformly + 0.2 mas for all model fits (see text). 
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ATLAS 12, and ATLAS 9 models. As a reference of the strength 
of the limb-darkening, UD (/ = 1 for < fx < 1 and 1 = oth- 



erwise) and FDD (/ = fi) model visibility functions are shown, 
with diameters 0ud and #fdd corresponding to our favorite 
PHOENIX model fit. Our measurements differ significantly from 
UD and FDD models, confirming the limb-darkening effect. 
All considered PHOENIX and ATLAS model CLV predictions 
lead to very similar model visibility functions up to the 2nd 
lobe, which are all consistent with our data. Hence, our data 
confirm the model-predicted strength of the limb-darkening ef- 
fect. Our measured values in the second lobe of the visibility 
function seem to lie systematically above the model predic- 
tions by ~ 0.5- lcr. It is not yet clear if these small differences 
are caused by systematic effects of our data calibration or by 
the model structures. The squared visibility amplitudes derived 
from the siderostat data show a systematic offset of about 1 <x 
towards larger values with respect to our best fitting curves. 
This is likely caused by possible small systematic calibration 
effects, since the calibration of our high squared visibility am- 
plitudes in the range ~ 0.8-0.95 derived from the siderostat data 
was difficult (see Sect. 14. 3> . 

The best fitting Rosseland angular diameters are well con- 
strained by our measurement. The measurements in the 2nd 
lobe of the visibility function also constrain the positions of 
the 1st minimum and 2nd maximum of the visibility function, 
which is a constraint of the diameter. This constraint is inde- 
pendent of possible small systematic calibration uncertainties 
of the \V\ 2 values. The reduced^ values for the different con- 
sidered models do not show significant differences, which is 
consistent with the very similar predicted shapes of the model 
visibility functions. These similarities are expected since our 
measurement is dominated by continuum photons , and the 



continuum forming region or the atmosphere is almost com- 
pact (see Fig. |3 . 

6. Discussion and conclusions 

Spherical PHOENIX models We have constructed a spherical 
hydrostatic PHOENIX model atmosphere for if/ Phe in Sect. 13.21 
Here, we confront this model's prediction for the CLV by com- 
paring it with our VLTI/VINCI measurement of the visibil- 
ity function in the second lobe. We find that the model pre- 
dicted shape of the visibility function is consistent with our 
VLTI/VINCI measurements. Simultaneously, the Rosseland 
angular diameter derived from the model and spectrophotome- 
try ff^*' = 8.0 + 0.4 mas agrees well with the Rosseland angu- 
lar diameter derived from the same model and our VLTI/VINCI 
measurements #r™ ci = 8.13 + 0.2 mas. These findings increase 
confidence in theoretical atmosphere modeling of cool giant (M 
III) stars. 

Spherical PHOENIX models with varied model parameters 
(r e ff, \ogg, M) that are still consistent with the values derived 
from spectrophotometry in Sect. 13.21 lead to the same best fit- 
ting values for 0r oss and xl- Hence, these model parameters 
cannot be further constrained by the shape of our measured 
visibility function up to the 2nd lobe beyond the constraints 
provided by the available spectrophotometry. 

The corresponding $ld and #r oss values for these differ- 
ent spherical models illustrate that #ld depends on the detailed 
model structure and that 6r oss is better suited to characterize 
the stellar diameter. The best fitting 6>ld values for the different 
considered spherical PHOENIX models differ by up to 0.24 mas 
or ~ 3%, while the 6fc oss values for the same models differ by 
only up to 0.02 mas or ~ 0.2%. 

Plane-parallel PHOENIX and ATLAS models Our plane- 
parallel PHOENIX and ATLAS models with the same model pa- 
rameters as used for the spherical model lead, as expected for 
a continuum-dominated measurement, to a very similar shape 
of the visibility function. These small differences can not be 
detected by our VLTI/VINCI measurements. 

The obtained Rosseland angular diameters derived by the 
plane-parallel models, which we assume to equal Sld, are con- 
sistent within our error-bars with the result obtained by our fa- 
vorite spherical model. However, there are systematic differ- 
ences of the obtained diameter values of up to ~ 0.1 mas or 
~ 1.5% as compared to the result obtained by our favourite 
spherical model. These differences can be explained by the 
different model geometries, line lists, opacity sampling tech- 
niques, and spectral resolutions of the employed models. 

Final parameter values for t^Phe Because of the 6% exten- 
sion of Phe's atmosphere, we consider our favorite spherical 
PHOENIX model the most reliable one. However, this can cur- 
rently not be verified by our VLTI/VINCI measurements. The 
diameter 9r oss = 8.13 ± 0.2 mas derived from this model and 
our interferometric data is the tightest available constraint on 
\fi Phe's diameter. From this angular diameter and our bolomet- 
ric flux from Sect. 12.21 an effective temperature of T e s- 3472 ± 



125 K is derived. 1ms means that the tightest constraint on 
\fi Phe's effective temperature comes from the model compar- 
ison with the available spectrophotometry in Sect. 13.21 which 
is r eff = 3550 + 50 K. With these values (6» R0 s S = 8.13 ±0.2 mas, 
r e ff= 3550 + 50 K) and the Hipparcos parallax, we derive a 
linear Rosseland radius of R - 86 + 3 R Q and a luminosity of 
logL/L =3.O2 ± 0.06. Together with the evolutionary tracks 
by Girardi et al. (2000 1, see Sect. l2.3l & Fig.^ these values are 
consistent with a mass of M — 1 .3 ± 0.2 M and a surface grav- 
ity of log g = 0.68^q}j. These values are summarized in the last 
Col. of Tabled They are consistent with the values derived by 
the different methods mentioned in earlier Sects. It would be 
an interesting further test to confirm the derived surface gravity 
and mass by means of a high-resolution spectrum. 

Future measurements Our VLTI/VINCI measurements 
could only probe the LD intensity profile of i^Phe averaged 
over the broad VINCI sensitivity band, and confirm the model- 
predicted strength of the limb-darkening. We have shown in 
Fig. [5] that the intensity profile and the strength of the limb- 
darkening effect of \fi Phe are expected to vary over our instru- 
mental bandpass, especially in narrow molecular bands. While 
we take this predicted variation into account for the computa- 
tion of the synthetic broad-band squared visibility amplitudes, 
it would be a better test of the model atmospheres to obtain 
spectrally resolved observations. Promising observations to 
better constrain atmospheres of cool giants would likely be di- 
rect limb-darkening measurements, i.e. observations with more 
than one resolution element across the stellar disk as used 
here, but at a number of well-defined narrow molecular and 
continuum bands with high spectral resolution. Theoretical 
studies seem to be needed and are planned to derive which 
observations in terms of wavelength bands and fundamental 
parameters of the target stars are best suited to constrain model 
parameters such as the model geometry and the treatment of 
molecular opacities. The required spectrally resolved limb- 
darkening observations can in the future be obtained with the 
upcoming scientific VLTI instruments AMBER and MIDI. 
AMBER will allow us to probe the stellar intensity profile with 
a spatial resolution A/B of up to ~ 1 mas and with a spectral 
resolution A/AA of up to 10000. AMBER can combine the 
light from three telescopes simultaneously, and hence obtain 
closure phases and triple amplitudes. The use of closure phases 
will likely enable the detection of additional surface features 
such as spots (Wittkowski et al. 2002). The use of the 1.8 m 
diameter ATs and the 8.2 m diameter UTs equipped with adap- 
tive optics will allow us to obtain these measurements with a 
high signal-to-noise ratio. 

7. Summary 

We have obtained VLTI/VINCI limb-darkening measurements 
of (/(Phe by probing the visibility function in the 2nd lobe. 
Our data confirm the strength of the limb-darkening effect as 
well as the Rosseland angular diameter as predicted by our fa- 
vorite spherical PHOENIX model atmosphere, the parameters 
for which were constrained by comparison to available spec- 



trophotometry and theoretical stellar evolution tracks. 1ms in- 
creases confidence in stellar atmosphere modeling of cool giant 
(M III) stars. 

We have derived high-precision fundamental parameters 
for tfrPhe, as listed in the last Col. of Tabled 

The measurements presented here are also a precursor for 
planned more detailed spectrally resolved limb-darkening mea- 
surements of a wider range of stars aiming at further constrain- 
ing the effects of model geometry and the treatment of atomic 
and molecular lines. Observations with the completed VLTI 
and its scientific instruments will enable these advanced studies 
of stellar atmospheres. 

Acknowledgements. Observations with the VLTI are only made pos- 
sible through the efforts of the whole VLTI team. We thank especially 
the VLTI commissioning team including the Telescope Instrument 
Operators for the operation of the interferometer. We would like to 
thank the entire PHOENIX model atmosphere team, in particular Peter 
Hauschildt and France Allard for helpful discussions. We thank Robert 
L. Kurucz for making his model atmosphere data publicly available, 
for providing us with a new ATLAS 12 model for our source i/fPhe, 
and for helpful discussions. We thank T. Szeifert for valuable com- 
ments and discussions on this project. We are grateful for the valuable 
comments by the referee M. Scholz which helped to improve this ar- 
ticle. This research has made use of the SIMBAD database, operated 
at CDS, Strasbourg, France. 

References 

Aufdenberg, J. P., Hauschildt, P. H., Baron, E., et al. 2002, Apl 570, 
344 

Aufdenberg, J. P., & Hauschildt, P. H. 2003, Proc. SPIE 4838, 193 
Ballester, P., Chavan, A. M., Glindemann, A., et al. 2002, Proc. SPIE 
4844, 300 

Baschek, B., Scholz, M., & Wehrse, R. 1991, A&A 246, 374 
Batten, A. H., Fletcher, I. M., & MacCarthy, D. G. 1989, Publications 

of the Dominion Astrophysical Observatory, 17, 1 
di Benedetto, G. P., & Foy, R. 1986, A&A 166, 204 
Bessell, M. S., & Brett, I. M 1988, PASP 100, 1134 
Bessell, M. S., Brett, I. M., Scholz, M., Wood, P. R. 1989, A&AS 77, 

1 (erratum A&AS 87, 621) 
Bessell, M. S., Brett, I. M., Scholz, M., Wood, P. R. 1991, A&AS 89, 

335 

Bevington, P. R., & Robinson, D. K. 1992, Data Reduction and Error 
Analysis for the Physical Sciences (WCB/MGraw-Hill), ISBN 0- 
07-911243-9 

Borde, P. I., Coude du Foresto, V, Chagnon, G., Perrin, G. 2002, A&A 
393, 183 

Burnashev, V. I. 1985, Abastumanskaya Astrofiz. Obs., Byull., 59, 83 
Burns, D., Baldwin, I. E., Boysen, R. C, et al. 1997, MNRAS 290, 
Lll 

Cardelli, I. A., Clayton, G. C, & Mathis I. S. 1989, Apl 345, 245 
Claret, A. 2000, A&A 363, 1081 
Claret, A. 2003, A&A 401, 657 

Cohen, M., Walker, R. G., Carter, B., et al. 1999, Al 1 17, 1864 
Coude du Foresto, V, Ridgway, S., Mariotti, I.-M. 1997, A&AS 121, 
379 

Coude du Foresto, V, Perrin, G., Ruilier, C, et al. 1998, Proc. SPIE 
3350, 856 

Coude du Foresto, V. 1998, ASP Conf. Series, 152, 309 

Crampton, D., & Evans, L. T. 1973, MNRAS 162, 11 

Davis, J., Tango, W. I., & Booth, A. I. 2000, MNRAS 318, 387 



Dyck, H. M., Benson, I. A., van Belle, G. T, Ridgway, S. T 1996, Al 

111, 1705 
Eggen, O. I. 1973, Apl 184, 793 

Evans, D. S., Menzies, A., & Stoy, R. H. 1957, MNRAS 1 17, 534 
Farge, M., 1992, Annual Review of Fluid Mechanics 24, 395 
Feast, M. W., Whitelock, P. A., & Carter, B. S., 1990, MNRAS 247, 
227 

Gezari, D. Y, Pitts, P. S., & Schmitz, M. 1999 Catalog of Infrared 
Observations, Edition 5 (Gezari+ 1999), VizieR On-line Data 
Catalog: 11/225. 

Girardi, L., Bressan, A., Bertelli, G., Chiosi, C. 2000, A&AS 141, 371 

Glindemann, A., et al. 2003, Proc. SPIE 4838, 89 

Grevesse, N., & Noels, A., 1993, in Origin and Evolution of 

the Elements, eds. N. Pranto, E. Vangioni-Flam & M. Casse, 

Cambridge: Cambridge Univ. Press, 14 
Hajian, A. R., Armstrong, I. T, Hummel, C. A., et al. 1998, ApJ 496, 

484 

Hanbury Brown, R., Davis, I., Lake, R. J. W., Thompson, R.I. 1974, 

MNRAS 167, 475 
Haniff, C. A., Scholz, M., & Tuthill, P. G. 1995, MNRAS 276, 640 
Hauschildt, P. H, & Baron, E., 1999, 1. Comp. Appl. Math 109, 41 
Hauschildt, P. H, Allard, F, Ferguson, I., Baron, E., Alexander, D. R. 

1999, Apl 525, 871 
Hofmann, K.-H, & Scholz, M. 1998, A&A 335, 637 
Houk, N. 1978, Michigan catalogue of two-dimensional spectral types 

for the HD stars, Vol. 2 (Ann Arbor: Department of Astronomy, 

University of Michigan) 
lankov, S., Vakili, F, Dominiciano de Souza, A., lanot-Pacheco, E. 

2001., A&A 377, 721 
lohnson, H. I. 1965, Comm. Lunar Plan. Lab. 3, 73 
Johnson, H. L., & Mitchell, R. I. 1975, RMxAA 1, 299 
Johnson, H. L„ & Mitchell, R. I. 1977, RMxAA 2, 269 
Jones, D. H. P. 1972, ApJ 178, 467 
Jones, D. H. P., & Fisher, J. L. 1984, A&AS 56, 449 
Kervella, P., Gitton, P., Segransan, D., et al. 2003, Proc. SPIE 4838, 

858 

Kervella, P., Segransan, D., & Coude du Foresto, V. 2003, in prepara- 
tion 

Kurucz, R. 1993, Limbdarkening for 2 km/s grid (No. 13): [+0.0] to 
[-5.0]. Kurucz CD-ROM No. 17. Cambridge, Mass.: Smithsonian 
Astrophysical Observatory, 1993 

Kurucz, R. 1996, ASP Conf. Ser. 108, 160 

Kurucz, R. 2003, private communication 

Leinert, C, Graser, U., Waters, L., et al. 2003, Proc. SPIE 4838, 893 
Lunt, J. 1919, ApJ 50, 161 
Moro, D., & Munari, U. 2000, A&AS 147, 361 
Percheron, I., Richichi, A., & Wittkowski, M. 2003, Proc. SPIE 4838, 
1424 

Perrin G., Coude du Foresto, V, Ridgway, S. T, et al. 1999, A&A 345, 
221 

Perryman, M. A. C, ESA 1997, The HIPPARCOS and TYCHO cat- 
alogues, ESA SP Series vol. 1200, Noordwijk, Netherlands: ESA 
Publications Division 

Petrov, R. G., Malbet, F, Weigelt, G, et al. 2003, Proc. SPIE 4838, 
924 

Plez, B., Brett, J.M., & Nordlund, A. 1992, A&A 256, 551 
Quirrenbach, A., Mozurkewich, D., Buscher, D. F, Hummel, C. A., 

Armstrong, J.T. 1996, A&A 312, 160 
Richichi, A., Lisi, F. 1990, A&A 230, 355 

Schlegel, D. J., Finkbeiner, D. P., & Davis M. 1998, ApJ 500, 525 

Scholz, M. 1985, A&A 145, 251 

Scholz, M. 2001, MNRAS 321, 347 

Scholz, M„ & Takeda, Y. 1987, A&A 186, 200 



Scholz, M., & Tsuji, T. 1984, A&A 130, 1 1 

Segransan, D., Forveille, T., Millan-Gabet, C. P. R., Traub, W. A. 

1999, ASP Conf. Ser. 194, 290 
Tango, W. J., & Davis, J. 2002, MNRAS 333, 642 
von der Liihe, O. 1997, in "Science with the VLT Interferometer" 

(Springer: Berlin), p. 303 
Wilson, R. E. 1953, General Catalogue of Stellar Radial Velocities 

(Washington, D.C.: Carnegie Institution) 
Wittkowski, M., Hummel, C. A., Johnston, K. J., et al. 2001, A&A 

377, 981 

Wittkowski, M., Scholler, M., Hubrig, S., Posselt, B., von der Liihe, 
O. 2002, Astron. Nachr. 323, 241 



